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Dipoles  in  Dissipative  Media 
by 

Ronold  W.  P.  King 

Gordon  McKay  Laboratory,  Harvard  University 
Cambridge,  Massachusetts 


Abstract 

The  general  problem  of  antennas  in  dissipative  media  is  reviewed  with 
reference  to  possible  applications  and  types  of  antennas.  The  question  of  bare 
and  insulated  antennas  is  discussed  with  reference  to  the  highly -conducting 
center-driven  cylinder  of  small  cross  section.  The  delta -function  generator  as 
a  convenient  idealization  of  feeding  by  transmission  lines  is  considered;  the 
apparent  difficulty  of  terminals  short-circuited  by  the  conducting  medium  is 
resolved. 

The  problem  of  determining  the  admittance,  the  distribution  of  current, 
and  the  elect romagretic  field  for  a  cylindrical  antenna  in  a  conducting  dielectric 
is  presented  with  emphasis  on  the  dual  requirements  of  reasonable  accuracy 
and  simplicity.  Methods  used  to  solve  the  integral  equation  for  the  current  in 
antennas  in  air  are  reviewed,  and  the  limitations  of  their  results  for  the 
calculation  of  the  complete  electromagnetic  field  are  pointed  out. 

An  approximate  method  of  solving  the  integral  equation  for  the  current 
in  an  antenna  in  a  dissipative  medium  is  described,  and  a  simple,  reasonably 
accurate  solution  is  displayed.  The  specific  evaluation  of  currents  and 
impedances  for  half-wave  dipoles  and  electrically  short  ar.tennas  is  summarized. 
The  separation  of  the  current  into  a  principal  part  for  vd.ich  the  electromagnetic 
field  may  be  evaluated  in  closed  form,  and  a  part  that  may  be  neglected  in 
important  special  cases  is  explained. 
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The  complete  electromagnetic  field  of  ?.  half-wave  dipole  in  a  dissipative 
medium  is  determined  first  in  cylindrical  coordinates,  then  in  spheroidal 
coordinates.  The  nature  of  the  wave  fronts  and  the  polarization  of  the  electric 
field  are  considered  fairly  close  to  the  antenna  and  at  greater  distances.  Similar 
results  for  the  electrically  short  antenna  are  also  given-  The  power  dissipated 
in  finite  parts  of  the  medium  is  determined. 

Introduction 

The  single  antenna  and  arrays  of  antennas  in  a  dissipative  medium  are 
aralytically  and  experimentally  interesting  and  offer  difficult  problems  that  have 
numerous  applications.  These  include,  for  example,  communicating  from  a 
submerged  submarine  or  from  a  mine  shaft,  and  carrying  oat  geophysical 
explorations  with  antennas  directly  imbedded  in  the  earth  or  lowered  into  air- 
fillea  holes.  Of  more  recent  interest  is  the  use  of  antennas  as  probes  in  the 
study  of  ionized  regions  and  plasma  sheaths  by  rr^eans  of  rockets.  In  general, 
the  regions  involved  do  not  consist  of  a  single,  homogeneous  and  isotropic 
medium-  There  may  be  two  distinct  layers  such  as  the  sea  oi  the  earth  and  the 
air  above  it;  alternatively,  as  m  the  interior  of  the  earth,  or  in  a  plasma  sheath 
in  the  ionosphere,  the  medium  may  have  continuously  varying  properties  or  it 
may  be  stratified  in  one  way  or  another.  Experimental  and  theoretical  studies 
of  bare  and  insulated  antennas  and  arrays  in  and  over  dissipative  and  dielectric 
media  that  are  isotropic,  have  continuously  varying  properties,  or  are  strati¬ 
fied  in  a  regular  manner  are  in  progress.  However,  this  xeport  is  limited  to 
a  study  of  the  properties  of  single  antennas  of  moderate  length  in  an  infinite, 
homogeneous,  isotropic,  dissipative  medium. 

Investigations  that  relate  to  antennas  in  dissipative  media  may  be 
separated  into  two  groups:  (a)  those  which  are  concerned  only  with  the 
electromagnetic  field  at  large  distances  from  a  source,  so  that  in  the  interest 
of  simplicity,  the  antenna  may  be  idealized  far  beyond  practical  availability; 

(b)  those  that  treat  the  circuit  properties  of  the  antenna  itself  as  well  as  the 
field  that  it  generates.  The  source  used  in  the  former  group  usually  consists 
of  a  Hertzian  dipole  or  infinitesimal  doublet  that  is  represented  mathematically 
by  a  periodically  varying  electric  moment  concentrated  at  a  point. 
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Physically,  such  a  doublet  may  be  visualized  as  made  up  of  equal  positive  and 
negative  charges  oscillating  back  and  forth  in  opposite  directions  along  a  line 
in  the  limit  as  the  charges  become  infinite  and  the  length  of  their  path  vanishes. 
In  the  second  group  are  studies  of  thin  cylindrical  antennas  and  of  biconical 
structures  immersed  in  dissipative  media.  The  antenna  may  be  in  direct 
contact  with  the  conducting  medium  or  it  may  be  surrounded  by  an  intervening 
layer  of  insulation.  Note  that  a  sufficiently  short  antenna  ie  the  physically 
realizable  approximation  of  an  infinitesimal  doublet,  insofar  as  the  distant 
field  is  concerned. 

An  early  study  of  the  radiation  of  a  Hertzian  dipole  immersed  in  a 
dissipative  medium  was  made  by  C.  T.  Tai  [1],  who  obtained  the  electro¬ 
magnetic  field  and  the  Poynting  vector  in  the  well-known  forms  for  doublets  in 
air,  but  with  the  permittivity  and  propagation  constants  both  complex  instead  of 
real.  He  noted  that  the  total  power  transferred  into  a  dissipative  medium,  as 
obtained  from  an  integration  of  the  normal  component  of  the  Poynting  vector 
over  a  sphere  with  its  center  at  the  doublet,  becomes  infinite  when  the  radius 
of  the  sphere  is  reduced  to  zero  and  concluded  that  "it  is  impossible  to  speak 
of  the  total  power  radiated  by  a  Hertzian  dipole  when  the  latter  is  in  direct 
contact  with  a  dissipative  medium.  "  He  then  proceeded  to  analyze  the  doublet 
enclu&ed  in  an  insulating  sphere  [2],  Actually,  the  power  transferred  across  a 
spherical  surface  that  encloses  charges  oscillating  along  a  line  is  not  obviously 
defined  in  the  limit  as  the  radius  of  the  enclosing  sphere  vanishes.  Indeed,  it 
is  shown  in  a  later  section  that  the  power  radiated  by  an  electrically  short 
dipole  in  a  dissipative  medium  cannot  be  obtained  oy  integrating  the  normal 
component  of  the  Poynting  vector  over  a  spherical  surface.  It  may  be  remarked 
in  passing  that  the  infinitesimal  doublet  has  been  used  as  an  idealized  source  in 
numerous  fairly  recent  studies  of  the  electromagnetic  field  in  a  conducting 
half-space  [3,4,  5,  6],  since  the  properties  of  a  finite  radiating  syscem  could  in 
this  way  be  avoided. 

The  integral  equation  and  its  formal  solution  for  the  current  in  a 
cylindrical  antenna  of  finite  length  immersed  in  a  dissipative  medium  was 
formulated  by  Tai  [7]  in  a  manner  paralleling  the  analysis  of  King  and 
Middleton  [8]  for  an  antenna  in  the  air.  The  essential  difference  is  that 
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the  previously  real  permittivity  and  propagation  constant  have  become  complex, 
with  a  resulting  complication  of  the  kernel  and  of  integrals  that  occur  in  the 
iteration.  Owing  to  a  lack  of  tabulated  functions,  Tai  did  not  evaluate  his  formal 
solution.  Indeed,  the  analytical  difficulty  associated  with  the  infinite  admittance 
of  the  delta -function  generator- -which  was  not  well  understood  at  the  time-- 
led  him  to  drop  further  work  on  the  cylinder  and  turn  his  attention  to  the  more 
tractable  problem  of  the  insulated  biconical  antenna  in  a  dissipative  medium. 
Tai's  work  on  the  dipole  has  been  extended  somewhat  by  Macrakis  [9], 

Harrison  [10],  and  Harrison  and  Denton  [11]  who  made  approximate  evaluations 
of  the  impedance;  more  recently  King  and  Harrison  [12],  and  King,  Harrison, 
and  Denton[13],  have  carried  out  complete  analyses  of  the  circuit  properties 
respectively  of  the  half-wave  dipole  and  of  the  electrically  short  antenna  in  a 
dissipative  medium.  These  studies  are  based  on  the  approximate  method 
proposed  by  King  [14]  for  the  solution  of  the  integral  equation  for  the  current  in 
a  cylindrical  antenna.  The  present  investigation  is  directed  to  the  determination 
of  the  circuit  properties  of  a  thin  cylindrical  antenna  of  moderate  but  arbitrary 
length  and  to  a  consideration  of  the  complete  electromagnetic  fields  generated 
by  the  currents  in  such  antennas. 

Review  of  the  Theory  of  Cylindrical  Antennas  in  Air 

Since  the  analysis  of  the  properties  of  an  antenna  in  a  dissipative 
medium  is  a  considerable  complication  of  the  problem  of  the  same  antenna  in  a 
perfect  dielectric,  it  is  well  to  review  briefly  the  extensive  theory  of  the 
cylindrical  antenna  in  air.  A  simple,  physically  realizable  circuit  consists  of 
a  cylindrical  conductor,  center-driven  from  a  balanced  two-wire  transmission 
line,  as  shown  in  Fig.  la.  Since  the  currents  in  the  antenna  and  in  the  line 
satisfy  two  simultaneous  integral  equations,  their  determination  is  a  formidable 
problem.  If  the  two  conductors  of  the  transmission  line  are  very  close  together, 
the  significant  interaction  of  the  line  and  the  antenna  is  confined  to  a  small 
region  near  their  actual  junction.  For  the  line,  this  may  be  approximated  by  a 
reactive  network  of  lumped  elements  characteristic  of  the  line  combined  with  the 
impedance  Zq  of  the  antenna  also  as  a  lumped  element,  as  shown  in  Fig.  lb. 
From  the  point  of  view  of  the  antenna,  the  driving  field  across  the  end  of  the 
feeding  line,  vdiich  is  distributed  over  a  short  length  of  the  antenna,  may  be 
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treated  as  an  impressed  field  that  is  concentrated  in  a  ring  around  the  center 
of  the  antenna,  in  the  form  E  (z)6(z)  ,  where  6  (z)  is  the  Dirac  delta  function. 
Clearly,  this  emf  iz  equivalent  to  a  discontinuity  in  scalar  potential  across  a 
pair  of  knife  edges  separated  ty  a  vanishingly  small  distance  at  z  =  0.  Such  an 
idealized  generator  at  the  center  of  the  antenna  evidently  includes  an  infinite 
capacitance  across  the  knife  edges,  so  that  the  input  susceptancc  must  also  be 
infinite.  It  was  shown  by  Wu  and  King  [15]  that  in  principle  the  infinite  current 
associated  with  this  capacitance  may  be  subtracted  out.  Moreover,  since  it  is 
confined  to  an  extremely  short  distance  adjacent  to  the  knife  edges,  it  is  in 
practice  automatically'  omitted  from  the  total  current  when  this  is  determined 
approximately  by  any  method  of  solution  that  represents  the  current  by  a  few 
terms  in  a  series  of  continuous  functions.  Thus,  the  practical  problem  may  be 
approximated  by  an  isolated  cylindrical  antenna  with  a  delta- function  generator 
at  its  center,  as  shown  in  Fig.  ic.  The  impedance  of  this  antenna,  after  the 
knife-edge  current  has  been  subtracted  out  or  omitted,  is  the  lumped  load  for  a 
transmission  line  with  a  suitable  terminal-zone  network.  The  nature  of  the 
lumped,  corrective  networks  required  for  different  connections  to  various  types 
of  lines  is  discussed  elsewhere  [16]. 

The  serious  attempt  to  determine  the  distribution  of  current  in  a  thin 
cylindrical  antenna  by  analytical  means  rather  than  by  assuming  it  empirically, 
as  is  still  commonly  done  in  the  so  called  emf  method,  begins  with  the  work 
of  L.  V.  King  [17]  and  especially  of  Hallen  [18]  whose  integral  equation  is  the 
basis  of  most  modern  theories.  For  a  perfectly  conducting  tube  of  very  small 
wall  thickness  and  radius  a  that  extends  from  z  =  -h  to  z  =  h  ,  the  integral 
equation  may  be  expressed  in  the  form 

h 

4t rp  1  A  (z)  =  \  I  (z1)  K  (z  ,  z')  dz' 
o  z  j 

-h 

=  lH2L[Ccoskoz  +  sinko|z|]  (1) 

bo 

where,  as  shown  in  Fig.  lc,  £(z’)  is  the  total  axial  current  at  z'  . 
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If  I{z)  includes  currents  on  the  inner  and  outer  surfaces  ox  the  tube,  it 
follows  that 

I(h}  -  C  (2) 

at  the  center  of  the  antenna  is  that  of  an  idealized  delta  function  generator. 
A^  (z)  is  the  vector  potential  at  the  surface  of  the  cylinder.  It  is  in  the  ix>rentz 
gauge  that  satisfies  the  Sommerfeid  radiation  condition.  The  kernel  K  is 
given  by 


K 


?  -jk  R 


,  n  1  C  e'^o 
(z’z)  =  2?  J  R 


,  -jk  R, 
d0  =  e  °  1 


R 


-IT 


where  R  =  z-z1)^  +  ( 


,  .  61.  2  B  f.  .,2  ,  2 

2a  sx ,  Rjs/iz-z-}  +a 


(3) 

(4) 


As  usual,  u  is  the  permeability:  h  =  I20ir  ohms,  if  is  the  characteristic 
impedance  of  free  space;  kQ  =  to/c  is  the  wave  number  C  is  a  constant  to 
be  determined  from  the  boundary  condition  (2). 


Approximate  solutions  of  integral  equations  substantially  like  (1)  have 

been  obtained  by  Hallen  and  others  by  a  method  of  iteration  which  depends 

implicitly  on  the  observation  that  the  ratio  A  (z)/I(z)  of  the  vector  pi  tential  to 

z 

the  current  along  the  antenna  is  approximately  constant  and  predominantly  real. 
Solutions  of  this  type  appear  in  the  form 


!(*) 


j2irV 


e 


r  v 

Jo 


M  (z)  4-  M,  (z}/v  +  M~(z)/y^  +  •  ’ 

Ox  Ct 

Fo(h)  +  Fjih)/  V+  F2(h)/  V2  +  •• 


(5) 


where 

Mo(z)  =  sinkQ(h  -  |z|)  ,  FQ(h)  =  cos  •  (6) 


The  first-order  terms,  M^(z)  and  Fj(h)  may  be  expressed  in  terms  of 
generalized  sine  and  cosine  integrals,  the  higher  order  terms  are  more 
complicated  and  must  be  evaluated  by  numerical  methods.  Unfortunately,  even 
with  the  most  sophisticated  definition  of  the  expansion  parameter  V  ,  at  least 
a  second-order  solution  is  required  in  order  to  obtain  quantitative  accuracy 
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for  antennas  with  electrical  lengths  in  the  range  0  <  k  h  <  27 r  .  A  solution  of 
this  type  ic  that  of  King  and  Middleton  [8. 19]  in  which  V  is  defined  to  be  the 
absolute  value  of  the  function 

h 

Y{z)  =  M^zMKIz.z’Jdz1  (6) 

-h 

at  the  point  where  M^fs)  has  its  greatest  value  in  the  range  0  <  z  <  h  . 

Extensive  ''omputations  of  the  impedances  [19]  of  cylindrical  antennas  in  the 
range  of  electrical  half-lengths  given  by  0  <  kQh  <  7  show  the  second-order 
King-Miadleton  values  to  be  in  good  agreement  with  experiment.  The  distri¬ 
butions  of  current  for  selected  lengths  have  also  been  computed,  but  the  second- 
order  formula  is  far  too  complicated  for  the  convenient  evaluation  of  electro¬ 
magnetic  fields.  For  this  purpose  the  rather  crude  zero-order  term  has  been 
used. 

The  quantitative  accuracy  of  second-order  results  calculated  'rom  (5)  in 

the  King-Middleton  ^rom  have  been  verified  theoretically  on  the  one  hand  by  the 

comparable  variational  solutions  of  the  integral  equation  (1)  by  Storer  [20]  and 

Tai  [21];  and,  on  the  other  hand,  by  the  detailed  study  of  the  integral  equation 

with  Fourier  series  methods  by  Duncan  and  Hinchey  [22].  These  investigators 

converted  the  integral  equation  into  a  set  of  simultaneous  equations  with  the 

Fourie-  coefficients  of  the  current  distribution  as  the  unknowns  and  then  carried 

out  calculations  to  the  25tn-order.  The  numerica.l  results  for  k  h  =  ir/2  and  tt  and 

o 

with  h/a  =  60  ar.d  500 ir  differ  by  only  about  2  percent  or  less  from  the  King- 
Middleton  values.  Significantly,  even  with  a  solution  of  such  high  order,  there 
is  still  no  sign  of  the  large  current  associated  with  the  knife  edges  at  the 
driving  point  of  the  delta- function  generator. 

The  reason  for  the  rather  slow  convergence  of  the  series  in  (5)  even 

with  a  rather  carefully  selected  expansion  parameter  has  been  traced  by  Wu  [23] 

to  the  method  of  evaluation  of  the  arbitrary  constant  C  in  terms  of  the  boundary 

condition  (2).  Unfortunately,  it  is  precisely  at  z  =  h  v/here  the  ratio  of  vector 

potential  to  cunent  departs  most  from  the  assumed  constant  value --it  actually 

becomes  infinite  at  this  point.  It  is  the  zero-order  form  of  C  ,  namely  , 

C-.---ft.ank  h.  which  leads  to  the  zero-order  distribution,  sink  (h  -  Izl)  . 

2  c  *  o'  1  1 


TR336 


-8- 


It  turns  out  that,  for  the  input  current  and  the  current  along  most  of  the  antenna, 

ve 

a  better  value  of  C  would  be:C=-w°  tank  (h  +  6}  where  6  is  a  small  length 
that  increases  with  the  radius  of  the  antenna.  Evidently,  such  a  value  of  C 
leads  to  a  small  nor.-vanishing  zerc-order  current  at  z  -  h  in  violation  of  (2). 
Primarily  owing  to  the  rather  poor  approximation  in  the  zero-order  value  of  C, 
the  solution  (5*  even  in  second-order  does  not  maintain  its  accuracy  as  the  length 
of  the  antenna  is  increased.  In  order  to  handj.e  specifically  the  very  long  an¬ 
tenna,  Wu  [23]  developed  an  asymptotic  solution  of  the  integral  equation  (1) 
bets^d  on  the  Wiener-Hopf  technique  and  a  method  of  evaluating  C  which 
properly  locates  the  distribution  of  current  along  the  antenna,  instead  of  re¬ 
quiring  it  to  vanish  at  the  end.  The  impedance  of  long  antennas  has  been 
computed  from  the  new  formula  by  Beaton  and  Wu  [24]  for  electrical  lengths  up 
to  k  h  =  30  .  The  new  theory  does  not  yield  a  simple  expression  for  the  current 
for  use  in  the  evaluation  of  electromagnetic  fields.  However,  the  radiation 
field  can  be  obtained  directly  from  the  Fourier  transform  of  the  current. 


The  major  characteristics  of  the  radiation  field  of  dipoles  m  air  may  be 
determined  with  reasonable  accuracy  from  the  simple  sinusoidal  distribution 
that  is  the  leading  term  in  the  iterated  solutions  for  the  current.  However,  this 
zero-order  current  has  the  serious  defect  that  its  value  at  the  driving  point  is 
at  best  a  rough  approximation  of  the  correct  input  current.  Moreover,  if  it  is 
used  to  calculate  the  radiated  power  with  the  Poynting-vector  theorem,  the 
result  may  be  in  error  by  as  much  as  50  percent,  if  it  is  assumed  to  apply  to 
center-driven  antennas  with  practically  significant  radii.  It  follows  that  even  if 
the  field  patterns  calculated  from  a  sinusoidally  distributed  current  are  an 
acceptable  approximation,  the  power  apparently  supplied  to  the  antenna  at  its 
terminals  and  the  power  radiated  may  be  grossly  in  error  and  mutually 
inconsistent.  If  the  input  admittance  is  known  accurately  from  a  higher-order 
theory,  these  difficulties  are  not  serious  for  antennas  in  air,  since  the  correct 
total  power  supplied  is  then  available  and  it  is  known  without  further  calculation 
that  this  is  equal  to  the  power  radiated.  When  an  antenna  is  immersed  in  a 
conducting  medium  the  problem  is  much  more  complicated,  since  power  is 
dissipated  throughout  the  medium. 
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What  i o  required  for  the  quantitatively  accurate  representation  of  the 
circuit  and  field  properties  of  dipole  antennas  and  arrays  is  a  formula  for  the 
current  that  is  sufficiently  simple  to  permit  the  evaluation  of  the  field  and  at 
the  same  time  sufficiently  accurate  to  yield  good  approximations  of  the  input 
admittance  and  of  the  radiated  power.  Such  a  formula  has  been  derived  by 
Storer  [10]  by  variational  methods  and  by  King  [14]  with  a  modified  iterative 
procedure.  Since  the  former  is  not  conveniently  applied  to  more  than  one 
antenna,  and  a  future  study  of  coupled  antennas  in  dissipative  media  is  contem¬ 
plated,  attention  is  focussed  on  the  latter,  which  is  as  useful  for  parallel  arrays 
as  for  a  single  antenna. 

The  required  relatively  simple  formula  for  the  current  is  obtained  as  the 
approximate  solution  of  a  rearranged  form  of  the  integral  equation  (1),  viz. , 

h 

4iru  lA  (z)  -  A  (hVl  =  \  I(z')K.(z.  z')dz' 

’  O  "  Z  *  ‘  z  ‘  J  G 

-h 


,i47T 


where  the  new  kernel  is 


Kd(z,z!)  =  K(z,  z')  -  K{h,  z ') 

and 


(7) 


(8) 


U 


Az(h)  = 


h 


I 


-h 


I(z')  K  (h,  z')dz  . 


The  shorthand  symbols 


and 


F  =  cos  k  z  -  cos  k  h- 
oz  o  o 


M  =  sin  k  (h  -  |z  [  ) 
oz  o'  II' 


(9) 


(10) 

(11) 


are  used,  together  with  Fq  (h)  which  is  defined  in  (6).  The  advantages  of  the 
rearranged  form  (7)  of  the  integral  equation  (1)  are  several.  The  integral  is, 
as  indicated,  proportional  to  the  vector  potential  difference  rather  than  to  the 


4 
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vector  potential  itself,  so  that  it  vanishes  at  z  =  h  as  does  the  current.  The 
right-hand  member  of  the  equation  is  the  sum  of  two  terns  that  are  individually 
related  to  well-known  distributions  of  current  and  vector  potential  difference. 

The  shifted  cosine,  F  ,  is  a  dose  approximation  of  the  current  and  the  vector 
potential  difference  along  an  unloaded  receiving  antenna  in  the  plane  wave  front 
of  a  distant  transmitting  antenna.  The  sine  term,  Mqz  ,  is  the  zero-order 
distribution  along  a  center-driven  antenna;  it  is  much  more  exactly  the  distribu¬ 
tion  of  both  current  and  vector  potential  along  an  ideal  two-wire  transmission 
line  with  an  open  end  at  z  =  h  and  a  delta-function  generator  at  z  =  0.  These 
facts  suggest  that  F  may  be  interpreted  as  an  approximation  of  the  distribu¬ 
tion  of  current  or  of  the  vector  potential  difference  that  is  maintained  by  the 
interaction  of  the  more  widely  separated  elements  in  an  antenna,  whereas  Mqz 
is  the  distribution  maintained  by  a  generator  when,  as  in  the  transmission  line, 
there  is  no  significant  interaction  between  widely  separated  current  elements. 


This  interpretation  is  confirmed  by  the  fact  that  in  (7)  'viQz 
coefficient  Ye  ,  the  actual  driving  voltage  of  the  generator,  whereas  the 
coefficient  U  of  Fqz  is  proportional  to  that  part  of  the  vector  potential  that 
has  been  subtracted  out  on  the  left  because  it  is  active  along  the  entire  antenna. 
In  general  and  as  a  first  approximation,  a  concentrated  generator  excites  a 


current  with  the  distribution  M 
distribution  F 


oz 


a  distributed  field  excites  a  current  with  the 


An  approximate  solution  of  (7)  may  be  obtained  if  the  integral  is  separ¬ 
ated  into  two  parts  of  which  the  one  varies  like  M  ,  the  ether  like  F 

r  oz  ’  oz 

This  separation  is  easily  accomplished  by  inspection,  once  the  kernel  has  been 
expanded  as  follows: 


Kd(z,  z')  =  KdR(z,  z')  +  jKdI(z,  z') 


«»  \ 
\lCai 


where 


KdR,(z>  z') 


7 r  | 

J_  c 

r cos  k  R 

cos  k  R, 

cosk  R, 

rifl1  -  °  1 

cosk  R-- 

1  o 

o  h 

o  in 

Ztt  J 

-rr 

L  R 

r\ 

nh  - 

d0  “  Rj 

Rlh 

(12b) 
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and 


V 

r 

sink  R 
o 

sinkoRh" 

d8'  =  - 

sink  R, 
o  1 

sink  R,, 
o  111 

j 

-71’ 1 

R 

u 

Rh  J 

L  Ri 

'  Rih  , 

(12c) 


The  subscript  h  on  R  and  Rj  denotes  the  values  defined  in  (4)  with  z  =  h  . 
Since  R  and  Rj  become  very  small  and  K^^(z,z')  correspondingly  very  large 
when  z'  is  near  z  ,  it  follows  that  the  principal  contributions  to  the  part  of 
the  integral  that  has  K^^(z,  z')  as  kernel  come  from  elements  of  current  very 
near  z’  =  z  .  This  means  that  the  part  A^„(z)  of  A^(z)  that  depends  on 
KdR(z,z')  varies  like  I(z)  .  On  the  other  hand,  since  K^j(z,z')  is  very  small 
when  z'  is  near  z  ,  the  principal  contributions  to  that  part  A  T(z)  of  A  (z) 
that  depends  on  K^j(z,  z')  come  from  all  the  elements  of  current  that  are  at 
some  distance  from  z  .  If  it  is  now  assumed  that  the  current  is  the  sum  of 
two  parts, 

I(z)  =  Ia(z)  +  iy(z)  (13) 

which  by  definition  have  the  leading  terms 

V>  ~  F0Z  i  Vz>  ~  Moz  <14> 


it  is  clear  that  (7)  may  be  separated  reasonably  into  the  following  parts: 


j  [lU{z,)Kd(a,»,)+ jly^'JK^a^'Mda' 


(h) 


02 


(15) 


I 

-h 


1  y  (  z'lKdR  (z>  z0  dz' 


. 

'  WV 


M 


02 


(16) 


Each  of  these  integral  equations  may  now  be  solved  by  iteration  in  the  King- 
Middleton  manner,  and  their  solutions  for  I^(z)  and  Iy(z)  combined  in  (18) 
to  give  I(z)  .  The  formula  for  the  current  so  obtained  may  then  be  substituted 
in  (9)  in  order  to  evaluate  the  constant  U  in  terms  of  Vq  .  The  result  is 


I(z)  = 


SoydR 


sink  (h  -  |z  !)  +  T(h)(cos  kQz 


coak  h) 
o  ' 


cos  k  h 
o 


+  higher  order  terms  (17) 
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The  indeterminate  form  obtained  when  kQh  =  7r/2  may  be  evaluated  to  give 
jZirVe 

I(z)  =  7  v/'  ~  [  sink  | z  |  -  1  +  T'  (z)  cos  k  z]  +  higher  order  terms  (18) 

VdR  o-o 

The  real  expansion  parameter  and  the  complex  constant  T(h)  are 

expressed  in  terms  of  tabulated  sine,  cosine,  and  exponential  integrals.  They 
are  functions  of  h/a  and  kQh  .  Explicit  formulas  and  numerical  values  are  in 
the  literature  [14,  25].  The  higher  order  terms  in  (17)  and  (18)  involve  essentially 
the  same  integrals  as  the  corresponding  terms  in  the  King -Middleton  solution 
(5),  but  their  contributions  to  I(zf  ,  at  least  in  the  first-  and  second-order 
terms,  are  very  much  smaller.  It  follows  that  the  quasi-zero-order  terms  in 
(17)  and  (18)  have  an  accuracy  that  lies  somewhere  between  the  first-  and  the 
second-order  solution  in  the  form  (5)  in  the  range  0  <  kQh  <  3ir/2  .  Moreover, 
the  error  is  largely  in  the  susceptive  part  of  the  current  in  a  small  region  very 
near  the  driving  point.  As  a  specific  example,  consider  a  half-wave  antenna 
with  h/a  -lb  .  The  distribution  of  current  as  obtained  from  (18)  is  shown  in 
Fig.  2  .  The  corresponding  admittance  and  impedance  are 

[Y  ]  =  (9.  07  -  j4.  67)  x  10  ^mhos  :  [Z  ]  -  82.  8  +  i39.  2  ohms  . 

The  second-order  theoretical  values  (which  are  in  good  agreement  with 
experiment)  are: 

[Y  ],  =  (9.  38  -  j4.  52)  x  10  ^mhos  :  [Z  ]7  =  86.  5  -!-  j41 . 7  ohms  . 

The  values  obtained  by  the  emf  method  independent  of  h/a  are 

Y  .  =  (10.  22  -  j5.  94)  x  10  ^mhos:  Z  .  =  73. 13  +  j42.  5ohms  . 
emf  J  ’  emx  J 

As  compared  with  the  second-order  value,  the  admittance  given  by  (18)  for  the 
half-wave  dipole  is  about  5  percent  in  error  in  the  conductance,  3  percent  in  the 
susceptance:  the  corresponding  errors  in  the  values  obtained  by  the  emf  method 
are  9  percent  and  31  percer.c.  For  greater  lengths  the  results  of  the  emf 
method  deteriorate  very  rapidly. 

It  may  be  concluded  that  (17)  and  (18)  combine  reasonable  accuracy  with 
simplicity  in  both  the  distribution  of  current  and  the  admittance. 


FIG.  2  QUASI -ZEROTH -ORDER  CURRENT 
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The  Cylindrical  Antenna  in  a  Dissipative  Medium 
i.  The  General  Case 

The  physical  problem  under  consideration  is  to  determine  the  current 
in  a  highly  conducting  cylindrical  antenna  immersed  in  a  homogeneous, 
isotropic,  dissipative  medium  of  great  extent.  The  antenna  is  center-driven, 
for  example,  from  a  shielded-pair  line  that  lies  in  the  neutral  plane.  As  with 
the  antenna  in  air,  it  is  convenient  to  approximate  this  configuration  by  an  ideal 
line  with  the  impedance  of  the  antenna  and  a  suitable  terminal-zone  network  as 
a  lumped  load.  For  this  purpose  the  antenna  is  imagined  center-driven  by  a 
delta -function  generator.  In  this  case,  the  knife  edges  of  the  generator  term¬ 
inals  are  equivalent  to  an  infinite  admittance  that  includes  both  a  capacitance 
and  a  conductance.  The  infinite  current  that  is  associated  with  the  charging  of 
the  knife-edge  capacitance  and  that  crosses  from  one  edge  to  the  otiier  by  way  of 
the  dissipative  medium  may  again  be  subtracted  out  in  principle,  since  it  is 
confined  to  an  extremely  short  distance  on  each  side  of  the  generator.  In 
practice,  it  is  excluded  from  a  solution  that  approximates  the  current  in  the 
antenna  by  a  series  of  continuous  functions. 

The  integral  equation  for  a  perfectly  conducting  cylindrical  antenna 
immersed  in  an  infinite,  homogeneous,  isotropic,  dissipative  medium  and 
center-driven  by  a  delta- function  generator  is  formally  like  the  equation  for  the 

<te 

same  antenna  in  air  if  the  complex  dielectric  factor  §  =  -j  —  is  substituted 
for  eQ  and  p  replaces  pQ  .  It  is  assumed  that  p  is  real.  This  is  equiva¬ 
lent  to  the  replacement  of  the  real  wave  number  kQ  =  by  the  complex 

propagation  constant 

k  =  P  -  ja  -  w/pe^  /I  -  jp  =  w/pe^  lf(p)  -  jg(p))  (19) 

where  p  =  <r  /we  and  f(p)  Jh(i  sinh  *p)  ,  g(p)  =  sinh(isinh  *p)  . 

Tables  of  the  functions  f(p)  a*.d  g(p)  are  available  in  the  literature  [26].  In 
addition,  the  real  characteristic  impedance  f  =  J\x  /e  =  wp  /k  is  replaced 
by  the  complex  value 

r  -  /5  -  “fi  -  f 
4  /?  k  Sae 


(20a) 
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where 

5a  •  -£*=  ,  i- =*»“  f  •  (201 

In  these  relations  and  cr^  are  the  real  effective  permittivity  and  con¬ 
ductivity,  respectively.  In  terms  of  the  complex  permittivity  e  =  e'  -  je'1 
and  the  complex  conductivity  a  =  o-1  -  jo-"  the  real,  effective  values  are  given 
by 

c  =  e  e  =  e'  ,  cr  =  c r'  +  (21 

e  o  er  u>  ’  e  1 

In  the  rearranged  form  (7)  the  integral  equation  is 


$  •'■ggjg  [ukrkz  ♦  KMks] 


Uk  =  Sf  j'i(d')Kt(h>z')d,' 

-h 

F,  (h)  =  cos  kh  -  cos  p  h  cosh  a  h  +  j  sin  p  h  sinh  a  h 
k 

F.  -  coskz  -  cos  kh 
kz 

=  (cospzcoshaz  -  cosph  cosh  ah)  +  j(sinpz  sinh  az  -  sin  P  h  sinh  ah) 

M.  =  sink(h  -  |z  |) 


sinP(h  -  |z  j)  cosh  a  (h  -  |z  | )  -  j  cos  P  (h  -  |z  | )  sinh  a(h  -  jz  |) 


The  kernels  are  given  by 


Kkd(z,z')  =  Kk(z,z')  -  Ky.lh,*')  -= 


-ikR,  -jkR,, 
,,  .  e  “  1  e  J  lh 


The  difference  kernel  may  be  separated  into  two  parts  as  follows: 
KkR(z,z')  =  y*,*')  +  jKkI(z,z') 


(29a) 
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where 

cos  p  R.  coshaR,  cos  |3  R,  coshaRu 

KkR<z> 2’>  - - i: - 1 - t - *  <28i> 

1  h 

K^z,  z1)  =  -R^  [  sinP  R^(coshaR,  -  sinh  aR^ )  +  j  cos  f3R^  sinh  aRj] 

+  Rj^ [  sinpRj^(cosh  qR^  -  sinhaRj^)  +  j  cospR^sinh  a  R^]  (21c) 

Note  that  Kkj(z,z')  is  not  real  when  a  /  0;  it  does  reduce  to  the  real 
K^j(z,z')  whe  a  =  0  .  This  equation  (22)  has  a  much  more  complicated 
kernel  than  (7)  and  the  right  member  is  also  more  involved,  since  the 
distribution  terms  F^z  and  Mkz  are  complex,  and  contain  the  additional 
parameter  a  . 

The  type  of  solution  desired  for  the  current  in  an  antenna  in  a  dissipative 
medium  is  one  corresponding  to  (17)  for  the  antenna  in  air  that  provides  good 
approximations  of  both  the  distribution  of  current  and  the  admittance  in  a  form 
that  is  sufficiently  simple  to  permit  the  direct  integration  of  the  integrals  for 
the  electromagnetic  field.  It  is  in  order  to  obtain  such  a  solution  that  the  form 
(22)  of  the  integral  equation  was  chosen.  The  approximate  solution  of  this 
equation  may  be  carried  out  in  a  manner  closely  paralleling  the  procedure 
described  for  the  antenna  in  air.  Indeed,  the  kernel  (27)  has  already  been 
separated  into  two  parts  such  that  the  one  Kkp(z,z')  ,  is  very  large  when  ?. 1 
is  near  z  so  that  the  principal  contributions  to  Az^(z)  are  from  the  currents 
near  z  ;  and  the  other,  K^j(z,  z1)  ,  is  relatively  small  near  z'  =  z  so  that  the 
principal  contr;buiions  to  Azj(z)  are  from  currents  at  some,  distance  from  z  . 
As  before  in  (13),  the  current  may  be  expressed  as  the  sum  of  two  terms  and 
the  integral  equation  appropriately  separated  into  two  parts  as  in  (15)  and  (16): 

C  j47TkUk 

J[iu(2')Kkd(z»s,)  +  %(z,)Kki<z>z,)ldz'  ^i^pyh)  Fkz  <29> 


r  j2irkV 

J  Iv<*,>KkR<z’»,>dB,‘i*^  M: 


(30) 
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Thess  equations  may  be  solved  separately  by  iteration  and  the  resulting  solutions 

0 

added  to  obtain  1(a)  in  terms  of  Vq  and  may  then  be  expressed  in 

terms  of  with  the  substitution  of  I(z)  in  (24)  .  The  solution  for  the  current 
is 


j2irkVe  [higher 

I(z)  = - — - .  °  .-v-  [  sink(h  - \z  j)  +  T.  (h)  (cos  kz  -  cos  kh)] +<  order  (31) 

up^kpcoskh  '  (terms 

where  is  a  complex  expansion  parameter  and  T^(h)  a  complex  constant 

that  depends  on  h/a ,  ph  ,  and  a/p  .  Unfortunately,  the  integral  functions 
that  occur  in  T^(h)  are  not  in  general  available  in  tabulated  form.  However, 
explicit  evaluation  has  been  carried  out  for  the  half-wave  dipole  [12]  in  a 
moderately  conducting  medium  such  that  ph  =  tt/2  and  ah  <  1  and  for  the 
electrically  short  dipole  [13]  in  an  arbitrary  medium  defined  by  ph  <  0.  3  , 
a  <  P  . 

2.  The  Half- Wave  Dipole 

For  the  half-wave  dipole,  the  approximations  coshah  =  1  and 
sinhah  =  ah  are  made  and  (31)  reduced  to  the  form 


Ji5l_  =  zj.27r(L.-  > 


V°Se~ 
o  ■  er 


^  °  [sin  P  |  z  [  -  1  +  T^  ( ^)  cos  pz]  +  higher  order  terms  (32) 

mVkR  K 


where 


.  [33) 

The  normalized  current  I(z)/ve  Jk  is  shown  in  Fig.  3  with  I(z)  =  I"  i  jl ' 

o  e  r  z  z 

as  a  function  of  z  with  2a/p  =  j  /tie  ■£  as  parameter.  The  normalized 
_  e  o  er _ 

impedance  Z  and  admittance  Y/  J eer  are  shown  in  Fig.  4.  Both  the 

distributions  of  current  and  the  impedance  are  in  agreement  with  experimental 

results  obtained  by  Iizuka.  It  may  be  added  that  measurements  have  been  made 

with  larger  values  of  c  /fcJeQc  than  those  used  in  the  theoretical  calculations 

represented  in  Fig.  4.  These  indicate  that  the  resistance  curve  in  Fig.  4 

reaches  a  maximum  as  tr  /we  e  is  increased  further  and  then  bends  down  to 

e'  o  er 

approach  zero.  The  corresponding  distributions  of  current  change  from  concave 

outward  as  in  Fig.  3  to  concave  inward  as  tr  /u>e  e  is  increased. 

°  c  o  er 


-5 


l'z/V  V«7 


0 


10  x  I0'3 
Amp/Volt 


I  \/V4Tt 


FIG.  3  NORMALIZED  CURRENT  I2  =  Iz+jli  IN  HALF- 
WAVE  DIPOLE  IN  MEDIUM  WITH  CONDUCTIVITY 
c  AND  DIELECTRIC  CONSTANT  c  =  €0er 


FIG.  4  IMPEDANCE  AND  ADMITTANCE  OF  A  DIPOLE  WITH  ELECTRICAL 
HALF  LENGTH  /3h=-rr/2  WHEN  IMMERSED  IN  A  MEDIUM 
WITH  CONDUCTIVITY  cr  AND  DIELECTRIC  CONSTANT  e  =  <-0cr  • 

/5  =  a;  'f/17  =  2ir/  X 
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3.  The  Electrically  Short  Antenna 

The  current  in  the  electrically  short  antenna  defined  by  ah  <  1 ,  (3h  <  1 
is  readily  obtained  from  (31)  as  a  series  in  powers  of  {3h  and  ah  .  The  first 
few  terms  are: 

«■>  = ^{(>  *  *¥)(■  -  ¥)  *  13fa2  -  ilA3  (>  •  £ 

The  corresponding  admittance  is 

Yo(k)  =  Go(k)  +  jB0(k)  (35a) 


where 

°«<k|  1  jj  h * ^Vr/3)  (‘  -  $  +3T^4i  I*  -  + 5“4)  j  <3i*> 

»  n.1  —JhJ*.  I.  a2).  §VfI,  .  4  q2  f  of)  .  p4h4  ,a/;  _  ,<,<£  afjl  ,35 

p2/'  3  r  v  v/  3<n-3,fi\  vv/r 


In  these  formulas 


F=1+%2fl  =  1+.1.08 


n-  3 


n-  3 


(36) 


and 


Also 


=  2ln(  iih/a )  . 


(37) 


(38) 


These  expressions  are  general  and  apply  to  all  values  of  a  and  (3  that  satisfy 
the  relations  ph  <  0.  3  ,  ah  <  0.  3  .  The  fact  that  it  is  necessary  to  retain 
terms  up  to  and  including  fourth  powers  of  the  small  quantity  ph  is  an  indica¬ 
tion  of  the  complicated  nature  of  the  admittance  of  an  antenna  in  a  dissipative 
medium  with  no  restriction  on  a  .  The  reason  is  obvious:  if  a  =  0  as  in  a 
perfect  dielectric,  the  leading  terms  are 


4  4 

r  .  3'rrp  h 
o'V  -  3?eVdR(  Si-  3) 


B  (k  ) 
o'  o' 


27T(3h 

V'dR 


(39) 
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On  the  other  hand,  when  a  =  p  as  in  salt  water, 


G  (k)  = 

°  e  cR 


B0<k) 


4  7rF63h3 

VdR 


(40) 


where  F  is  defined  in  (36) .  Note  that  GQ(k.)  in  (39)  is  a  pure  radiation  con¬ 
ductance;  it  is  very  small  compared  with  the  suscspcance.  In  (40)  GQ(k)  is 
determined  entirely  by  dissipation  in  the  medium,  the  contribution  from 
radiation  is  negligible;  GQ(k)  is  very  large  compared  with  BQ(k)  .  It  is  clear 
from  the  general  expression  (35a)  that  the  radiation  conductance  predominates 
for  only  a  very  small  range  of  a  near  zero.  This  indicates  that  the  short 
dipole  behaves  as  a  radiating  antenna  only  when  a  is  very  small;  for  larger 

values  it  is  essentially  a  pair  of  electrodes. 

Z  (k)  -<Q(k)  X  (k) 

The  normalized  impedance  ^ - +  j  — -  and  admittance 

Y  /WM?  =  n  +  iR  fkl*5  of  7k  short  anfptina  with  e  Bh  =  0.  3  and  Q  =  10 

o'" "e  “o'*  '  Je  *  o'*  ' ■’e 

are  shown  in  Fig.  5  as  a  function  of  a/p.  The  curves  for  RQ(k)/fe  and 
GQ(k)?o  are  also  shown  on  a  logarithmic  scale  in  Fig.  6  in  order  to  show  their 
extremely  rapid  rise  in  the  range  of  very  small  values  of  a/p  .  Note  that 
RQ(k)/?e  has  a  maximum  between  a/p  =  0.  5  and  0.  6  . 

If  the  medium  in  which  the  admittance  YQ(k)  is  measured  is  an  ionized 
region  such  as  the  ionosphere,  the  conductance  GQ(k)  and  the  susceptance 
B0(k>  may  be  related  directly  to  ,.ie  concentration  of  electrons  and  to  their 
collision  frequency.  Details  are  given  elsewhere  [13]. 

The  Electromagnetic  Field  of  an  Antenna  in  a  Dissipative  Medium 

1.  The  General  Case  and  the  Field  at  Distant  Points 

The  electromagnetic  field  of  a  cylindrical  antenna  of  half-length  h  with  a 
total  axial  current  I(z)  when  immersed  in  a  dissipative  medium  may  be  deter¬ 
mined  from  the  vector  potential 

r  -jkR 

Az(r)  =4¥  JI(z,)“R“dz'  (41) 

-h 


where 
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when  expressed  in  the  cylindrical  coordinates  p,  6  ,  z  .  The  formula  (41)  may 
be  used  to  determine  the  vector  potential  at  any  point  except  within  di 
the  ends  of  the  antenna  comparable  with  its  radius  a  .  At  distances 
from  the  center  of  the  antenna  that  are  large  compared  with  the  half-length  h  , 
that  is,  when 


2 

r 


R  =  r  -  z'  cos  9  , 


(43) 


an  approximate  formula  for  the  vector  potential  in  the  radiation  zone  may  be 
used.  This  is 


.  r.  .  u  e  C  n  jkz’cos  9 

Az(r )=-&?— T~  \  Kz’)eJ 


The  corresponding  non-vanishing  cylindrical  components  of  the  magnetic  and 
electric  fields  are  given  by 


8Az(r)  2  9  A  (r)  k2  dA^r) 

B0'r)  =  "  9p  »  EP  •r)  =  u  dpdjZ  ;  Ez(r)  =  "V"  p  8/a  (P  8p  )  •  <45) 


In  the  radiation  zone  the  significant  components  are  conveniently  expressed  in 
the  spherical  coordinates,  r  ,  0  ,  $  .  They  are  given  by 

B*(r)  =  -jkAz(r)  sin  0  ;  Eq  (r)  =  -juA^  (r)  sin  9  .  (4i 


In  order  to  determine  the  complete  electromagnetic  field  it  is  necessary 
to  evaluate  the  integral  in  (41)  or  at  least  the  derivatives  of  this  integral 
required  in  (45).  For  the  radiation  field  at  distances  that  satisfy  (43)  it  is 
sufficient  to  evaluate  the  simpler  integral  in  (44).  A  simple  and  reasonably 
accurate  representation  of  the  current  in  a  cylindrical  antenna  of  zero  to  mod¬ 
erate  length  is  that  given  in  (31).  Unfortunately,  if  this  distribution  is 


i&intm 
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substituted  in  (41)  and  (45)  is  evaluated,  one  of  the  three  terms  involved  cannot 
be  obtained  in  closed  form.  On  the  other  hand,  the  integration  of  (44)  is 
straightforward.  With  (46)  the  result  is 


7ee'jkr 

=  e;<7)  *  °  coVkh  [Fm(9,  kh)  t  Tk(h|Gm(6,kh)] 

kr 


(47) 


where 


i  kh 


F  _  cos(khcos  8)  -  cos  ] 

~  m'  *  '  ~  sin  0 

_  /a  ...  sinkhcos(khcos6)  cos  0  -  cos  kh  sin(khcosO) 
™'W»Kn'=  sin  0  cos  0 


m 


/aa\ 
\  / 


(49) 


These  expressions  can  be  separated  into  real  and  imaginary  parts  and  the 
electric  and  magnetic  fields  evaluated.  However,  since  the  amplitude  decrease 
as  e  /r  ,  it  is  evident  that  the  range  of  r  in  which  the  field  is  significant 
and  in  which  (43)  is  satisfied  is  not  great  unless  a  is  sufficiently  small  to 
satisfy  the  conditions 


2  2 

ah  «  1  ,  sinh  ah  =  ah  ,  cosh  ah  4  1  . 


(50) 


Nevertheless,  the  general  directional  properties  of  the  field  not  very  close  to 
the  antenna  may  be  determined  from  (48)  and  (49)  even  when  a  =  g  .  When  (50) 
is  satisfied  the  field  factors  defined  in  (48)  and  (49)  become  somewhat  simpler. 
Specifically, 


where 


Similarly, 


Fm(0,kh)  =  Fmr(0,kh)  +  jahFm.(0,kh) 


F  (0  kh)  cos(ghcos0)  -  cos  gh 
■t!mr'w,  '  -  sin0 

F  (0  kh)  -  8^n(P^cos  cos  Q  ~  singh 
nil'0’  '  "  sin  0 


Gm(0,kh)  =  Gmr(0,kh)  +  jahFmi(0,kh) 


(51a) 


(51b) 


(52a) 


i 
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where 


G  (0  kh)  -  s*n  Ph  cos  (fih  cos  0)  cos  Q  -  cos  gh  sin  { 
'nr'  ’  '  Sin  0  COS  0 


Gm.  (0,  kh)  =  sin  ph  sin  (Phc°s  0)  tan  0  . 


(52b) 


Also  in  (47), 


cos  kh  =  cos  gh  +  jah  sin  gh 


When  gh  =  ir/2  ,  (51)  and  (52)  become 


r  .  -7T 

sin  \-j  CO! 
=!  +  jah  )_ - 


in  (-5-  cos  Q)  cos  8-1 


Fm(0»iT  -  ja^  =  Fm(0 


G  (0,4tt  -  jah)  =  F  (0,-5)  +  iah  sin(-5cos  0)tan0 
m  ’  i.  m  ’  L'  -  i. 


■jj  here 


COS  (-yCOS  0) 

F  (0,4y)= - =-r -  (54c 

nv  ’  2 '  sin  0  v 

Since  from  (53)  coskh  =  iah  ,  it  follows  with  the  readily  verified  approximation 
(that  is  valid  only  when  gh  =  7r/2) 


tan0  sin(-^cosO)  = 


and  the  definition  (33),  that  (47)  becomes 


V" 


ve  e-(a  +  jg)r 

o _ 

y.  r 
kr 


where 


cos  0  -  sin(-ycos  0) 

Hm<Q>7>  = - slnTcosG -  =  <!-?)  Fm<0’T>  *  ^ 

IT 

This  is  the  complete  field  generated  by  the  current  (32).  Note  that  H  (6,-y) 
is  the  part  of  the  field  contributed  by  the  component  of  current  sing|z|-J  . 
The  factor  (1  -  jalnr/2)/  (1  -  jah)  in  (57)  as  compared  with  (32)  is  a  conse¬ 
quence  of  the  attenuation  of  the  field  over  the  small  distances  z'cosQ  in 
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R  =  r  -  z'cosQ  .  The  field  factors  F  (0,  —  )  and  H  (0,-r-)  are  shown  in 

m'  2 '  m'  5  2  ' 

Fig.  7, 


If  the  approximate  relation  on  the  right  in  {56b)  is  substituted  in  (56a), 
the  approximate  electric  field  is  given  by 


ve^"(a  +  JP)r 
Er  -  _I°f _ 

E0  "  V.  r 

kr 


1 


1  -  jah7r/2\ 

1  - jah  )  _ 


F  (0,  -y ) 
mv  ’  2 


This  field  would  be  maintained  accurately  by  the  current 


^j_=-j2^(l-j./p)|-1-^t  cospa> 

^  ^o  kr  L  J 


vevT~ 

ov  er 


(57) 


(58) 


and  only  approximately  by  (32).  Note  that  in  (32)  the  input  current  is  propor¬ 
tional  to  [-1  +  T^(^)]  whereas  in  (58)  it  is  proportional  to  [1  -^+TJ(^)]. 
Since  the  former  gives  the  correct  value,  the  latter  is  evidently  in  error.  The 
explanation  for  this  discrepancy  is  clear.  The  correct  current  in  (32)  yields 
the  correct  field  in  (56a)  and  the  correct  radiated  power.  If  the  current  is 
approximated  by  a  purely  cosinusoidal  distribution,  its  amplitude  can  be  adjusted 
so  that  the  input  current  is  correct  with  the  factor  [-1  +  T£.(^)]  ,  or  so  that 
the  field  and  the  power  radiated  are  correct  with  the  factor  [l  -  -jy  r  T^(^)  ]  . 

In  the  former  case  the  power  radiated  is  too  small,  in  the  latter  case  the  input 
current  is  too  large. 


Graphs  of  Hie  eiecliie  field  defined  in  (47)  when  ph  ~  IT,  (1  =  0,  h/a  -  75 
are  shown  in  Fig.  8.  It  may  be  concluded  from  the  forms  of  (54a)  and  (54b)  that 
the  field  patterns  when  ah  is  small  do  not  differ  greatly  from  those  when  a  =  0. 
Significant  differences  may  be  anticipated  when  ah  is  not  small  compared  with 
one,  in  particular,  when  a  is  of  the  same  order  of  magnitude  as  (3  .  For 
their  evaluation  the  complex  parameter  T^(h)  must  be  computed. 

2.  The  Field  Near  the  Dipole;  Approximate  Currents  and  Cylindrical  Coordinates 


In  order  to  determine  the  electromagnetic  field  nearer  the  antenna  than 
is  permitted  by  the  condition  (43)  it  is  necessary  to  carry  out  extensive  numer¬ 
ical  computations  for  individual  cases  or  to  make  approximations  in  the  distri¬ 
bution  of  current  in  order  to  obtain  an  integrable  form.  Fortunately,  the  latter 
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procedure  is  straightforward,  since  the  term  sink(h  -  Jzj)  in  the  genera] 
expression  for  the  current  (31)  is  integrable  when  substituted  in  (41)  and  (45). 
However,  it  is  not  to  be  expected  that  a  satisfactory  approximation  can  be 
achieved  if  the  entire  term  Tk(h)  (coskz  -  coskh)  in  (31)  is  simply  omitted. 

In  particular,  it  is  essential  that  the  total  power  supplied  -o  the  antenna  be  a 
good  approximation  even  if  the  distributin'-,  of  current  along  the  antenna  is  some¬ 
what  in  error.  This  may  be  achieved  by  a  reairangement  of  (31)  in  the 
following  manner: 

I(z)  =  [l(0  )/A{  0  )]  [A(z>  +  B(z)]  (59a 


where 


A(7.\  =  Fi  +  T  th\  f---~  £2^\  1 
'  •  L  K'  '  \  sinicn  /  J  L  coskn 


t<»)  [■ 


cos  kz  -  1  + 


/ 1  -  cos  kh 
l  sinkh 


k.|z|"|  . 


Note  that 


(59c) 


K  0  )  =  v“Yo(k)  = 


j2-jrkV 


A(0) 


where  YQ(k)  is  the  admittance.  The  two  parts  of  the  current,  A(z)  and  B(z) , 
are  chosen  so  that  B(  0  )  -  0  and,  hence,  contributes  only  to  the  distribution 
along  the  antenna  and  not  at  all  to  the  input  current  and  the  admittance,  which 
are  dec  rmined  entirely  by  A(  0  )  .  The  plan  is  to  neglect  B(z)  in  the  integra¬ 
tion  to  determine  the  electromagnetic  field  since  it  contains  the.  term  that  has 
not  been  integrated  in  closed  form.  The  anticipated  error  is,  therefore,  that 
made  by  the  omission  of  B(z)  .  Investigation  shows  that  A(z)  alone  is  quite  a 
good  approximation  when  (3h  <  37r/4  and  very  poor  when  ph  ~  ir  .  It  is 
particularly  good  when  (3h  <  tt/Z  and  the  maximum  current  is  at  the  driving 
point.  Indeed,  when  ph  =  tt/Z  and  (ah)  «1  ,  B(z)  =  T^  ( ^-)[cospz  +  sinP  |z  | -1]. 
The  contribution  by  this  component  of  current,  shown  shaded  in  Fig.  9  when 
a  =  0  ,  is  seen  to  be  quite  small.  Similarly,  when  ph  <  1  , 


,22  2 

B(z)  =  Tk(h)^-  (l-^)-{ 

h 


>-¥>• 


This  is  the  verv  small  difference 
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between  triangular  and  parabolic  distributions  when  required  to  have  the  same 
value  at  z  =  0  .  If  B(z)  is  neglected  in  (59a)  the  approximate  distribution  is 


simply 


1(a)  =  1(0)  ■Slrt-k3(^|z|)  ;  ph<3ir/4,  n<S 


where  1(0)  is  given  in  (59d) . 

When  (60)  is  substituted  in  (41)  the  vector  potential  itself  cannot  be 
obtained  in  closed  form.  However,  the  derivatives  required  in  (45)  may  be 
expressed  as  follows  [?. 7]: 

Be(7)=iHM  [e-JkRlh  +e-jl*2h  -  2coskhc'->kr] 

Ep(t)  =  e"jkRlh  +  *±Jl  e-jkR2h  .  ^coskhe-^] 

r  4"Kf  L  Rlh  R2h  r  J 

_  -iu)Ml(O)  f  e  j^Rlh  e  ■'kR2h  2  -jkr 

E2(r)  =  -Lfejr*  [~V—  4  -R^“  -  7coskhe 

where 


Rlh  =  /{Z  -  h)2  +  pZ  ,  R2h  =  /{z  +  h)2  +  P2  ,  r  =/z2+p2 


For  the  part  of  the  current  represented  by  (60)  these  expressions  are  accurate 
and  give  the  complete  electromagnetic  field  within  distances  of  the  cylindrical 
antenna  comparable  with  the  radius  a  .  Since  all  distances  are  measured  from 
the  axis  rather  than  the  surface  of  the  cylinder,  these  expressions  have  no  mean¬ 
ing  for  p  <  a  .  They  are  useful  particularly  for  determining  the  near  field  of 
antennas  with  ph  <  v/Z  . 

If  (61)  -  (63)  are  multiplied  by  e,)Ut  ,  the  real  part  is  the  instantaneous 
electromagnetic  field.  However,  the  formulas  so  obtained  permit,  no  simple 
physical  interpretation  other  than  that  the  ends  and  center  of  the  antenna  may 
be  regarded  as  the  origins  of  complicated  outward-traveling  disturbances.  A 
somewhat  simpler  picture  may  be  obtained  for  the  half-wave  dipole  in  spheroidal 
coordinates  and  for  the  electrically  short  antenna  in  spherical  coordinates. 


Milliamperes  per  Volt 


FIG.  9  CURRENT  DISTRIBUTION  ALONG  A  HALF¬ 
WAVE  DIPOLE.  THE  BROKEN  CURVE  IS 
I'0c°s  k0z. 


TR336 


-25- 


3.  The  Half-Wave  Dipole;  Spheroidal  Coordinates* 

V/hen  (3h  =  ir/2  the  normalized  distribution  of  current  in  (60)  is 

}\z)  _  cos  flz  cosha(h  -  jz  | )  .  sin  (3  |z  |  sinh  a(h  -  | z  j  > 

1(0)  “  coshah  ^  coshah 

The  associated  distribution  of  charge  per  unit  length  is  obtained  from 

?.'z)  =  Req(z)  +  jlmq(z) 


(65) 


(66a) 


where  in  normalized  form, 


p  c  0  (z)  _  (3  cos  (3z  sinh  a(h  -  |  z  j )  -  a  sin  ft  j  z  j  cosh  a(h  -  |  z  j ) 

L  1(0)  j  ”  p  cosh  ah 

t  I  c  3  (z)~j  _  sin  P  1  z  1  cosh  a(h  -  1  z  1 )  +  a  cos  |3  z  sinh  a(h  -  j  z  | )  | 
L  1(0)  J  "  "  i_  P coshah  J 


(66b) 

(66c) 


The  real  and  imaginary  parts  of  l(z)/l(0)  and  c  q(z)/l(0)  are  shown  in  Fig.  10 
for  nh  =  0  f  0.  1  0.  5  t  .  0  ,  and  ir/ 2  ;  c  is  the  velocity  of  light. 

Let  the  spheroidal  coordinates  k^  ,  k_  ,  and  $  be  introduced  where 
kg  =  a  /h  is  the  reciprocal  of  the  eccentricity  of  a  spheroid  with  semi-major 
axis  a^  and  with  its  foci  at  the  ends  of  the  antenna,  k^=  a.^/h  is  the  corre¬ 
sponding  quantity  for  an  orthogonal  hyperboloid  of  two  sheets,  and  <J>  is  the 
azimuth.*!  angle.  These  coordinates  are  related  to  the  cylindrical  coordinates 
by 


fi  -  hy^kg  -  1)  (1  -  k£)  z  =  bk^k,  ; 

they  are  illustrated  in  Fig  11.  Note  that  1  =  kg  =  co  ,  -1=  =  1  .  The 

spheroidal  components  of  the  electromagnetic  field  that  correspond  to  the 
cylindrical  components  (61)  -  (63)  are  quite  complicated  when  a  /  0  .  It  is 
convenient  to  express  the  instantaneous  value  of  each  component  as  the  sum  of 
three  terms  as  follows: 

3  3  3 

B*(r,t)=  £  B#.(?,t)  ,  Ee(?,t)=  £  Ee.(r,t)  ,  E^?,  t)  =  £  Ehi(?,t)(67) 

i  =  1  i  =  1  i  =  1 


Much  of  the  work  in  this  section  was  carried  out  originally  by  Mr.  S.  T.  Yu. 
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B>..{p,t)=  -*A  . 


cos(7rk^/2)  sin(ut  -  tt^/Z) 

/(kf  -  i)(i  -  kh 


Eel(r’t}  =  '?nA 


L  v  '"e  *'h' 

cos(?rkv/2)  sin(ut  +  ^  -  7rk  /2)  *| 


E,  ,  (r  ,  t)  =  5  A, 

h  lv  ’  a  i 


°  '  L 

<*1  /  0”L-  /?  \  »f  JL 

bill  ^  <«  M  /  »  } fJ  I *  XN.^  f  J 


kh 


)K-d 


(68a) 


(68b) 


(68c) 


A]  =  cosh(_il')  exp  ; 


BvZ(r’t}  =  ; 


•sin(wk^/2)cos(u)t  -  irke/2) 
,/fk2  -  1)11  -  k?) 


’  e 


il ' 


Ee2(r,t)=-?„A 


a  2 


sin  (irk^/2)  cos  (wi  +  ^  -  irk^/Z) 

fir? 


K-Kw-K' 


Eh2(r,t)=  ?gA2 


cos(irk, /2)  sin(ut  +  -  Trk  /2) 


i-2  •- -  ••- 


k“-kh)(k-ii 


t/rvx 

a  - 

a2  ■  27Th 


/"7raku\  /  iroKe\ 

sinhlTrii  expr~2F/  ; 


TTak. 


cos 


B$3(r  ,t)  =  pA3 


E.3<r'*>"f.A3 


2p  /  \  2P 

(Ut  ~  I7r'/ke  +  kh  _  l) 


/(kZe  -  1)(1  -  k2) 

kecos  [ut  +  ^  -  ■£,7rv/k2  +  kh  ”0 
/(k2  +  k2-  i)(kf-k2)(1  -k2) 


(68d) 

(69a) 

(69b) 

(69c) 

(69d) 

(70a) 


(70b) 


FIG.  10  NORMALIZED  CURRENT  AND  CHARGE  IN  ANTENNA  IN  DISSIPATIVE 
MEDIUM.  (  APPROXIMATE  THEORY) 


P/z  h  0  01  0.2 


0.3  0.4 


FIG.  II  ELECTRIC  FIELD 
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Eu,(r,t)=  S„A, 


/  . — ; - - - ^ 

k^cos  [wt  +  /  -  ^-ir /k“  +  -  lj 

/t2  4.  V2  -  1  \  /l-2  -  l-2'.  /b-2  -  1  \ 


^e  *  kh  “  ')  >  e  kh'-  '-ke  “  1  > 


A3=1£h  sinh  llpi  exP  ( 


/  ira  /£ 


ke  +  14'1 


Note  that  ph  =  .ir/ 2  ,  /  =  tan"*(a/p)  ,  £q  =  Ufa  /  ,/a2  +  p2  . 

The  third  part  of  the  field  as  given  in  (70a)  -  (70d)  may*  be  expressed  in 
the  spherical  coordinates  r  ,  0,  and  $  in  the  following  equivalent  form: 


=  filial  .1  P.F.)  sinh  <I±)  evn(-nrl 

5'r  ’  Z TT  r  sin  0  °mn  ,  exp(  ar) 


(71a) 


tano) 

Ll93(r  >  *•)  =  2.ir 


siSif  Pr/  sinh(l|)  exp(-ar) 


r3  (r,t)=  0 


These  rather  elaborate  expressions  for  the  electromagnetic  field  reduce 
to  a  simple  and  readily  interpreted  form  when  the  attenuation  of  the  medium  is 
rather  small.  Specifically,  if  the  following  inequality  is  satisfied: 

a/p  <  ah  «  1  (72) 

the  entire  electromagnetic  field  of  the  half-wave  dipole  is  well  approximated  by 
the  first  terms  (68a)  -  (68c)  in  the  sums  in  (67)  with  /  =  0  and  with 


iimi 


SS:  exp(-aaj 


Except  for  the  addition  of  the  exponential  term  in  the  amplitude  (73),  the  field 
given  by  (68a)  -  (68c)  with  /  =  0  is  like  that  for  the  same  antenna  in  air.  The 
components  B^(r*,  t)  and  Eg(r,t)  are  in  phase  with  each  other  and  both  are  a 
quarter  period  out  of  phase  with  E^(r  ,  t)  .  The  surfaces  of  constant  phase  are 
the  spheroids  kg  =  constant  and  these  expand  outward  in  such  a  manner  that 
the  end  point  of  the  semi-major  axis  travels  with  the  velocity  c  .  The  electric 
field  at  any  point  is  elliptically  polarized  as  shown  in  Fig.  11.  The  components 
B^(r,t)  and  Ee(r,t)  decrease  in  amplitude  with  distance  from  the  antenna  in 
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a  manner  that  approaches  l/r  at  great  distances;  the  component  E,  (r\  t)  , 

z  ^ 

on  the  other  hand,  decreases  in  a  manner  that  approaches  l/r  .  The  entire 
field  is  exponentially  attenuated  in  amplitude  with  the  semi-major  axis  ag  =  hke 
as  the  variable.  The  instantaneous  value  of  the  outwardly  directed  component 
Sh(r/t)  of  the  Poynting  vector  associated  with  the  electromagnetic  field  that 
satisfies  (72)  is 


SeI  (0)  cos2  (7Tkh/2)  exp  (  -  TTak^/f>)  2 

S,  (r,t)  =  — 5— y  - s -  -  - —  3in  -  irk  /2) 

4'b 


(74) 


where  f  =  wp/(3  .  The  time-average  power  transferred  across  a  spheroidal 

surface  with  semi-maior  axis  a  =  hk  is 

e  e 

Te  =|l2(0)R®  e‘2aae  (75) 

G 

where  Rq  =  73.  1  ohms.  Since  terms  of  the  order  of  magnitude  ch  hairo  been 
neglected,  the  exponential  in  (75)  reduces  to  unity  when  ag  approaches  h  ,  and 
(75)  reduces  to  the  value  for  air.  The  difference 


T 


el 


-  Te2  ‘I  >2  (O'  Ro  t * 


"2aaei  _  -2aae2 


(76) 


where  agj  and  a^  arc  the  semi-major  axes  of  two  spheroids  with  ae,  >  agj, 
is  the  power  dissipated  in  heating  the  slightly  conducting  medium  in  the  volume 
between  the  two  surfaces.  As  pointed  out  in  conjunction  with  (57)  and  (58)  the 
correct  power  when  the  approximate  distribution  of  current  (60)  with  (3h  =  ir/2 
is  used  is  obtained  with  the  value  of  1(0)  given  by  (58)  instead  of  the  correct 
,alue  given  by  \3>Z).  Alternatively,  the  correct  power  is  also  obtained  if  1(0) 
is  obtained  from  (32)  and  the  correct  input  resistant  Rq  is  substituted  i:.  (76) 
for  the  radiation  resistance  R®  as  computed  with  the  approximate  current. 

When  a  is  not  sufficiently  small  to  satisfy  (72)  with  (3h  =  it / 2  the 
complete  expressions  (67)  with  (o8a)  -  (70d)  must  be  used  for  the  electromag¬ 
netic  field.  These  cannot  be  combined  to  give  a  simple  picture  of  outward 
traveling  surfaces  of  constant  phase  at  all  points  outside  the  half-wave  dipole. 

It  is  not  even  convenient  to  combine  the  first  and  second  parts  of  the  field  into 
a  spheroidal  wave  and  treat  the  third  part  in  the  forms  (71a)  -  (71c)  as  a 


spherical  wave  originating  at  the  center  since  components  of  both  the  second 
and  third  parts  become  infinite  along  the  z  axis  and  only  their  sum  is  finite. 
However  if  a  relatively  small  region  in  the  vicinity  of  the  antenna  is  excluded 
the  complete  field  can  be  approximated  by  spheroidal  waves  similar  to  those 
which  obtain  when  the  attenuation  constant  is  small.  In  effect,  the  spherical 
wave  originating  at  the  centei  of  the  antenna  must  be  approximated  by  a 
spheroidal  wave.  It  is  readily  verified  that  if  the  quantity 

d  --  1  -  /~(l  -  kl)/kZe  [1 

is  small  enough  so  that  it  may  be  neglected  completely,  the  third  part  of  the 
field  given  by  (70a)  -  (70d)  may  be  approximated  by 

cos  (tat  -  irk  /2) 

B  (7,t)  =  pA  - - - 3.—  (7 

*  /2  2 
Ak  -  l)(i  -  kT) 


.(r,t)  =  ?  A, 


TA, 


cos(wt  +  <f>  -  irke/2) 

M  -  kjj>  (i  •  kj  i 

kh  cos  (wt  +  <f>  -  TTk^/2 
ke'/lke-khKke-  ‘> 


a  -U2L 

A3”  Zirh 


IT  1rCl^ 

sinh  <^)  exp  (  - 


It  is  to  be  noted  that  over  a  considerable  range  near  the  z  axis  k^  is  very 
near  one  so  that  d  =  0  .  In  the  vicinity  of  the  equatorial  plane  where  k^  is 
very  small,  the  excluded  range  is  determined  by  the  magnitude  of  kfi  ,  it  is 
indicated  in  Fig.  12.  Specifically,  for  example,  d  <  0.  1  when 

(1  -  k£)Ae  <0.19  . 

With  the  approximate  expressions  (78a)  -  (78d)  and  (68a)  -  (69d)  the 
following  combinations  may  be  made: 
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Dzcos(wt  -  Tk  /2 ) 

B*2(r  ,t)  +  B  3(r  ,t)  ==  -  ^  =p  exp{-  -  ?p  ) 

5  1(0)  D,cos(ut  i-  4  -  Tk  /2)  ™ 


-)  (79b) 


Ehl(?,t)  +  Eh3(r,t)=^ 


£.1(0)  F,  cos  (ut  +  4  -  TTk&/2) 


/(^-k^(ke2-l) 


exp  (  -  -"Vo-" )  (79c) 


where 


7rak,  Trk.  ira 

D,  =  sinh  (-gg^)  sin(-T2)  -  sinh  («) 


7T  ak,  Trk,  k, 

Fj  =  cosh  (— 2fp)  sin  (“2")  +  k~  sinh^2p^ 


With  the  additional  notation, 


h  h 

—  cosh  ( — 2g —  )  cos  (  ^  ^ 


;h,  ,Tkh 


(80a) 


F2  =  sinh (-—i  )  cos  (-g- )  < 

the  complete  field  may  be  obtained  by  combining  (79a)  and  (79b)  respectively 
with  (68a)  and  (68b),  and  (79c)  with  (69c).  The  results  are 


g  /f  ti  .  . 

"  Znh 


D  sin(wt  +  A  -  7rk  /2) 


/k^-Dd-k^) 


r  1NA/  A  , 

=: —  fiXP("  ~zf~] 


(81a) 


?  1(0)  D  sin(ut  +  4  +  A  -  irke/2) 


/ikf-k^iu  -k^r 


exp(  -  -j^-)  (81b) 


?  1(C)  Fcos(ut  +  4  +V  -  irk  /2)  Tftke  . 

tt  i?  _  _2 - - - - exp  (  -  "a"- )  (8lc) 

27Th  r-* - =; - S -  lY 


1.0 


FIG.  12  REGION  EXCLUDED  BY  kg  »  (I -k„ );  ROTATION 
OF  PRINCIPAL  AXES  THROUGH  ANGLE  0' 
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lt:Vo  Y*  Q 

D  =  /d?  +  D~  ,  A  =  tan_1(D-/D.)  (82a) 

A.  Lt  C*  1 

and 

F  -  /F2  +  F2  ,  V  =  -tan“1(F2/F1)  .  (82b) 

The  formulas  (81a)  -  (81c)  for  the  spheroidal  components  of  the  field  of  a  half¬ 
wave  dipole  (ph  =  tt /2  )  in  a  dissipative  medium  without  restriction  on  the 
conductivity  are  substantially  like  the  formulas  for  the  same  antenna  in  a  medium 
that  is  only  slightly  conducting.  However,  the  formulas  for  the  unrestricted 
medium  are  not  useful  in  a  region  near  the  antenna  where  the  approximation 
d  =■  0  is  not  a  good  approximation.  In  addition,  they  are  more  complicated 
owing  to  the  appearance  of  the  different  phases  A  and  V  in  the  expressions 
for  the  mutually  perpendicular  components  Eg(r,t)  and  E,(r,t)  of  the 
electric  field. 


The  time-average  component  of  the  Poynting  vector  perpendicular  to  the 
spheroidal  surfaces  kg  =  constant  is 

£ 

—  T2/m  D  COSj^  exp(  -  TTQk  /p) 

Sh(r)=  22  - 2"  ,  -  - *  <83) 

8.  b  -  ^)^T.  k8)|k8  .  „ 

The  total  time-average  power  transferred  across  a  spheroidal  surface  for 

which  k  >2.2  is 
e  — 


T 

e 


?eI2(0)  exp(-7rake/p)  n  D 

-i~ 


4ir(l  1-  a2/p2) 


2  dkh 


K 


(84) 


where  =  Wfi/P  .  Owing  to  the  complicated  form  of  D  as  a  function  of  k^ 
the  integral  in  (84)  has  not  been  evaluated. 

The  general  formulas  (81a)  -  (81c)  indicate  that  the  electromagnetic 
field  of  a  half-wave  dipole  in  an  arbitrary  dissipative  medium  may  be  inter¬ 
preted  in  terms  of  expanding  spheroidal  surfaces  except  very  near  the  antenna 
where  more  complicated  conditions  exist.  These  surfaces  expand  with  foci 
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at  the  ends  of  the  antenna  and  with  an  outward  phase  velocity  along  the  z  axis 
that  is  given  by 

vp  =  u/P  (85) 

where  (3  is  the  real  part  of  the  complex  propagation  constant  k  as  defined  in 

(19).  However,  although  the  phase  of  B^r",  t)  ,  Ee(?,  t)  ,  or  E^i?,  t)  remains 

constant*  as  each  spheroid  moves  out  along  a  hyperbola  defined  by  a  particular 

value  of  k^  ,  the  actual  phases  at  different  points  along  such  a  spheroid  are  not 

the  same.  That  is,  the  expanding  spheroids  are  surfaces  of  constant  phase,  but 

the  constant  value  on  each  spheroid  is  a  function  of  k^  .  Thus,  respectively, 

for  B.(r’,  t),  E  (r,  t)  ,  and  E.  (r  ,  t)  the  phases,  (ut  +  a  -  irk  /Z) , 
i  G  n  G  ^ 

(u)t  +  i  +  a  ■  irke/2),  and  (ut  +  <f>  +V  -  vk^/Z)  remain  constant  as  t  and  kfc 
increase  *ogether.  At  each  point  on  a  given  spheroid  k  =  constant,  the  phase 
of  E^(r,t)  differs  from  that  of  B^fr^t)  by  the  constant  angle  /  =  tan  (a/p)  . 
The  phases  of  the  two  mutually  perpendicular  components  of  the  electric  field 
differ  by  the  angle  (a  -  V)  which  is  a  function  of  k,  . 

In  a  medium  that  is  dissipationless  or  has  very  low  attenuation,  the 
electric  field  is  elliptically  polarized  with  principal  axes  tangent  and  perpen¬ 
dicular  to  the  spheroidal  surfaces  of  constant  phase  as  shown  in  Fig.  12.  In  the 
dissipative  medium  with  a  unrestricted,  the  principal  axes  of  the  ellipses 
described  by  the  electric  vector  are  rotated  from  these  tangents  and  perpen¬ 
diculars  by  an  angle  8'  that  is  a  function  of  k^  .  The  polarization  ellipses 
for  the  electric  fieid  and  their  angle  of  inclination  with  respect  to  the  direction 
of  E^(r  ,  t)  may  be  determined  if  the  time  is  eliminated  between  the  expressions 
for  Efi(r,t)  and  E^(r,t)  as  given  in  (81b)  and  (81c).  Let  these  formulas 
be  expressed  as  follows: 


$  -*■  -*■ 

This  is  strictly  true  for  B,(r,t)  and  E  (r  ,  t) ,  since  k  does  not  occur  in 

the  phase  term  A  .  It  is  not  actually  true  for  E^(r  ,  t)  owing  to  the  occur¬ 
rence  of  kfi  in  the  term  (k.^/k^)  sinh  (TTa/2p)  in  Fj  which  contributes  to  the 
phase  in  V  .  However,  this  term  is  small,  except  quite  near  the  ends  of  the 
antenna,  so  that  the  general  picture  is  not  altered  significantly. 
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Eg(r,t)  =  Eg  sin  ( \i  +  A  ) 


Eh{r,t)  =  Ehco5{Q+  V  ) 


£aI-;0}  D  exp  (- 7rake/2(3) 


-  k?  )(1  -  kj) 


?aI(0)  Fexp(-7rak  Jl$) 
27rh  r~2  2  2 

/(ke2-k^)(kf -i) 


f2  =  (tot  +  ^  -  7rke/2)  .  (90 

If  (86)  and  (37)  are  solved  for  sin  £2  and  cos  £2  ,  and  these  quantities  are  then 
squared  and  added  to  equal  one,  the  following  equation  is  obtained: 


E2(?,  t)  2Eh(?,t)  Ee(?,t) 
E?  +  E2  + 


sin  (V  -  A  )  =  cos  (V  -  a)  (91) 


This  is  the  equation  oi  an  ellipse  that  has  its  principal  axes  rotated  with  respect 

to  the  directions  of  E,  and  E  .  Let  the  direction  of  EJ  be  rotated  from  E, 

n  e  n  ii 

by  an  angle  0'  as  shown  in  Fig.  12  .  The  direction  of  E^  is  then  rotated  by 
the  same  angle  with  respect  to  Efi  .  The  new  components  E^  and  E^  are 
related  to  the  old  ones  by  the  formulas 

Eh  =  E^ cos  0'  -  E^  sine*  (92a; 

Eg  =  E^sin  8'  +  E^cos  0'  .  (92b] 


If  the  angle  9'  is  chosen  so  that 


2E,  E 

tan 2  0'  = — 5 - %  sin (V  -  A  )  , 

E  -  E  £ 
e  h 


i 
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the  new  equation  is 


|E£(r,t)f  |‘E;(r,t) 

.  tt i  +  r"  ip' 


=  cos  (  V  -  a) 


where 


cos2  6' 


2  sin(V  -  A)  sinQ1  cos  91 

E,  E 
h  e 


sin2  8' 


M  \  2 

e:  i 


■  2a, 

sin  8' 


2  sin  { V  -  A)  sin  8 1  cos  8 1 

E.E 
h  e 


cos2  8' 


The  equation  (94)  is  that  of  an  ellipse  with  semi-principal  axes  E^cos  (V  -  A) 

and  E^cos(V-A)  •  Note  that  when  V  =  A  ,  as  in  a  dissipationless  medium, 

A1  =  0  ,  E.'  =  E.  E'  =  E  .  It  can  be  argued  from  symmetry  or  determined 

>  n  n  ;  e  e  *’  ' 

directly  from  the  formulas,  that  8'  =  0  when  =  0  ,  1  ;  it  follows  that  the 
electric  field  is  still  linearly  polarized  parallel  to  the  antenna  both  along  the 
z  axis  and  in  the  equatorial  plane  just  as  when  immersed  in  air. 

4.  The  Electrically  Short  Antenna 

As  pointed  out  following  (59a)  -  (59d),  the  approximate  distribution  of 
current  (60)  in  the  form 

I(z)  =  1(0)  (1  -  jz|  ,/h)  (96) 

for  the  electrically  short  antenna  ((3h  <  0.  3 ,  ah  <  O.^J)  _is  a  particularly  good 
approximation  since  the  omitted  term,  B(z)  =  T^(h)  [(1  -  ^)  -  (1  -  *^*-)  ] 

is  very  small.  The  input  current  is  to  be  determined  from  the  more  accurate 
formula  (34)  .  Since  (96)  is  a  special  case  of  (60)  with  kh  small  in  its  teal  one 
imaginary  parts,  it  follows  that  the  rigorous  expressions  (61)  -  (63)  for  the 
electromagnetic  field  of  the  distribution  (60)  apply.  They  may  be  converted  into 
a  more  common  but  also  slightly  more  restricted  form  by  a  series  expansion  in 
powers  of  the  small  quantity  kh  and  the  quantity  h/r  ,  which  must  also  be 
assumed  small  .  The  distance  r  =  r  z ^  is  measured  from  the  center  of 
the  dipole  to  the  point  where  the  field  is  calculated.  The  approximations  include 


F 
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Rj^  =  r  -  hcos  9  +  (h  sin0)  /2r  (9' 

R£k  =  r  +  h  cos  0  +  (h  sin  0)  /2r  (9' 

coskh  =  1  -  k2h2/2  .  (9' 

With  these  approximations  the  complete  electromagnetic  field  in  the  spherical 
coordinates  r  ,  0  ,  has  the  familiar  form 


.  inl(0)h 
4lT 

u  r  J 

(98a) 

.  jupKOjh 
47Tk 

(98b) 

.  up l(0)h 
4irk 

i5‘.3le"jkrsine  ■ 

(98c) 

Note  that  up/k  =  Jp/%  =  ?  .  In  this  derivation  terms  of  the  order  of  magnitude 

|kh|4  ai\d  h2/r^  and  higher  powers  have  been  neglected.  This  means  that  the 

field  given  in  (97a)  -  (97c)  is  not  valid  in  the  immediate  vicinity  of  the  antenna 

as  are  the  more  general  expressions  (61)  -  (63).  tf  desired,  the  above  formulas 

may  be  expressed  in  terms  of  the  equivalent  electric  moment  p  with  the 

z 

relation 


z)  dz  =  1(0)  h 


For  a  short  dipole  in  air,  k  =  k^  is  real;  wp/k  =  .  The  instantaneous 

field  is  obtained  from  the  real  parts  of  (93a)  -  (98c)  when  expressed  in  polar 
form  and  multiplied  by  e^Wt  .  Specifically, 


B<^r  ~  4t r 


poI(0)h  /l  +  k2  r2 


-  cos  (ut  -  k  r  +  tan  k  r)  sin  0 

2  o  o 


(100a) 


5  l(0)h  J\  -  kV  +  k4r4 


o _ o_ 

3 


cos[ut  -  kQr  +  tan  (kQr  -  l/kQr)]sin0  (100b) 
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^  ?  I{0)h  v4+k2r2 

Er(r  ,  t)  =  -2— — - ^ —  cos  t  -  k  r  +  cot  K  r)  cos  9  *  (lQOc) 

k  r 
o 

These  expressions  do  uoi  in  general  peimit  the  simple  interpretation  of  spher¬ 
ical  surfrees  of  constant  phase  that  travel  radially  outward  with  a  definite  phase 
velocity  as  is  true  in  terms  of  spheroidal  waves  of  the  half-wave  dipole  in  air. 
Each  component  may  be  so  interpreted,  but  the  phase  velocities  of  all  three 
components  are  different  until  k  r  becomes  sufficiently  great  in  the  radiation 

2  2  i  o  _  j 

zone  so  that  k  rL  »  1  .  tan_1k  r  =  ir/2  ,  cot  k  r  =  0  .  The  radiation-zone 
o  '  o'*  o 

field  is 


B^(r  ,t)  =  Eg(?,t)/c  = 

M0>koh 

47Tr 

sin  (u)t  -  kQ  r)  sin  0 

(10la) 

?I(0)k  h 

= 

O  O 

2irr2 

cos(ut  -  kor)cos  8=0  • 

(101b) 

These  formulas  represent  a  true  spherical  wave.  All  components  expand  with 
the  same  phase  velocity  c  ;  Ik,  (r  ,  t)  and  E„  {r  ,  t)  are  in  phase  with  each 

*  I*  — *•  —p  j*  — ► 

other  and  a  quarter  period  out  of  phase  with  E  (r,t)  .  B.(r  ,  t)  and  Efi(r  ,t) 
decrease  as  l/r;  E^(r,t)  decreases  as  l/r  so  that  when  r  is  large  this 
component  is  insignificant  as  indicated  in  (ldb)  . 

The  radial  component  of  the  complex  Poynting  vector  for  the  field  in 
(100a)  -  (100c)  is 


_  lK0)|2koh2?c 


S1,(r )  = 


32  v 


J. _ 1 _ 

?.  .35 

r  k  r 


2 

sin  8  , 


(102) 


and  the  total  time -average  power  radiated  is 
2  tt  T 


T  =  Re 
e 


j  d$>  Jsr(r)r2  sin9  d0  =  j  lZ  (0)  R® 


(103) 


where  the  radiation  resistance  is 
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?  k2h2  .  ? 

Re  -  — : —  =  20  h  ohms 
o  o-rr  o 


(1C4) 


Since  approximate  current  (96)  was  used  tc  determine  the  power  radiated 

0 

instead  of  the  mere  accurate  current  (34),  Rq  in  (104)  is  only  an  approximation 
of  the  actual  input  resistance  obtained  from  (35) .  The  more  accurate  value  has 
as  its  leading  term, 

?  k2  h2  V, 

R  (k  )  r  (1C5) 

o'  O  07T  ft  -  3  ' 


where  Vjr  =  2(£n  ^  -  1)  and  ft  =  2fn(2h/a).  When  ft  =  iO,  h/a  =  75,  the  ratio 
Vjr/(  ft-  3)  =  0.  95,  so  that  (104)  is  in  error  by  about  5  percent.  The  error 
decreases  as  h/a  increases. 

When  the  medium  is  dissipative  with  k  =  |3  -  ja  .  the  components  of  the 
field  may  be  expressed  as  follows: 


4>  '  4tt  r  r  2 

L  r 

jr  I?)  -  tJP1(°)h  j  (i  a _ SL 

9  47T  ji  \r  •  ((JZ  +  „ 


-ar  -i6r  .  , 

e  e  J  sin  I 


U  2  \ 

■P  ~  a  ) 

2  22  * 


2,  2  lo2  .  2.2  3 


2  .  2.  2 


(P"  +  a  )  r  '  (p  +  a  )r 


(p2  +  a2)2r3 


1  -ar  -jpr  .  , 

|  e  e  Jr  sine 


(10ba) 


(106b) 


„  uul(0)h  r 
Vr)  =  ~Ztt  ! 


P  ,  _ 2ag _  (  q 

,r2  2.  2  +  ..2  2.2  3  +JL2±  2, 

^(P  +  a  )r  (P  +  a  )  r  MP  +  a  ) 


,ft2  .  2,2  3  j 

(p  +  a  )  r  ' 


l2  -  a2  \  -ar  -jpr 

i  I  a*  o  i 


JK*  o 

e  e  Jr  cos0 


(106c) 


The  real  part  of  the  complex  Poynting  vector  is 


Re  Sr(r)  = 


1 1(0)  j2  p2h2  ? 
327T2  r2 


6  i *  ,  &  «  i  Zar 
-  (1  +  p  L)e 
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where 


L,  = 


2  .  ? 

(i  t  ..  /p“)pr 


WP 


„  2  /a2. 2  D2  2 
(1  +  a  /p  )  P  r 


?  7  ?  ?  -a 

(i  +  q-/p  rrr 


(108 


and  where  5  =  wp/p  .  Thf  .otal  time -average  power  transferred  outward 

across  a  spherical  surface  u.  radius  r  is  the  integral  of  (107)  over  the  surface 
of  the  sphere.  The  result  is 


T 

e 


lifo)|2e2h2f  _ 

IT? - +  pL)e 


(109) 


As  in  the  case  of  (103)  this  formula  is  an  approximation  that  is  quite  good  only 
for  extremely  thin  antennas.  It  may  be  corrected  to  apply  to  antennas  of 
somewhat  greater  thickness  by  multiplying  (109)  by  the  factor  Vjr/(  ft-  3). 

It  is  important  to  note  that  (109)  cannot  be  used  to  determine  the  radiation 
resistance  ol  the  snort  antenna  in  a  dissipative  medium  since  it  is  not  possible 
to  reduce  the  spherical  surface  across  which  the  power  is  transferred  to 
an  envelope  that  encloses  only  the  antenna  and  no  part  of  the  medium.  The 
smallest  spherical  surface  that  contains  the  antenna  has  the  radius  r  =  h.  but, 
since  in  the  transformation  from  cylindrical  to  spherical  coordinates,  terms  of 
the  order  of  magnitude  (h/r)  have  been  neglected  in  comparison  with  unity, 

(109)  is  meaningful  only  when  r  >  h  .  Since  a  sphere  of  radius  r  >  h  contains 
a  significant  part  of  the  dissipative  medium,  thp  power  transferred  across  its 
surface  is  not  the  total  power  radiated  from  the  dipole.  Note  in  particular,  that 
the  limit  r  =  C  is  meaningless. 

The  actual  power  supplied  to  the  antenna  may  be  obtained  from  the 
driving-point  admittance  Y  (k)  in  (34)  in  the  form  Tg  =  jV^G^k)  =  jl(0)|  Rj[k). 

Since  the  complete  expressions  for  G  (k)  and  BQ(k)  are  rather  long,  it  will 
serve  for  purposes  of  illustration  to  consider  the  two  special  cases,  a  «  p~ 
and  a  =  p  ,  for  which  the  formulas  are  much  simpler.  The  resistance  RQ(k) 
is  given  by 


Rc.(k)  = 


o  no 

'  o’  ' _ 

G2(k)  +  B2(k) 
o'  o 


(110) 
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When  a2«  p2,  G2(k)  «  B2(k)  so  that 

G  (k)  ?e(32h2V,  r 

Rc{k)  =  B2(k)"  &i'a:  3)  f 

O  *- 


r  +  6a(  ^-3) 

L  P2h 


(111) 


2  2 

if  only  the  leading  terms  are  retained.  When  a  =  p  ,  GJk)  »  B^(k)  ,  so  that 
from  (40) 


Vk>  5  G 


__  feV  JL 

(k)  “  47T  ph 


The  power  per  unit  input  current,  Py  ,  that  is  dissipated  within  a  sphere 
of  radius  r  >  h  is  defined  as  follows  in  terms  of  the  power  Te  supplied  to  the 
antenna  and  the  power  Tgr  transferred  outward  to  the  region  beyond  the  sphere: 


7JT  -  t  j 
‘~e  'er' 


p2  fc2  ?  yA 


r  |K0)| 

When  a2  «  p2  this  becomes 


ft  =  Ro<k>  • -RFTnTf  <1+?«e 


P  sAV  a  r, 

r  "  6tt(  u  -3)  6 


s.  rS>J£.-JJ.  -  Le'2arl 
P  ph  **  J 


114M 


where 


T  — 2 

L  "  Pr  +  „3  3 
p  r 


(114b) 


When  a  =  p  , 


P  _  ^dr 

r  47T 


-^(i  +  L)e 


(115a) 


where 


i  _  ,J_  .  _JL_  ,  i  , 

L  ~  (pr  +  2  2  +  3  3  }  ■ 

p  r  Zp  r 


(115b) 


It  is  instructive  to  consider  a  numerical  example.  Let  the  power 
dissipated  within  tl  e  radian  sphere  pr  =  I  be  determined  for  an  antenna  with 
ph  =  0.  3  and  h/a  =  75  .  When  a2  «  p2  ,  e  2ar  =  e  4  1  -  2a/p ;  L  =  4  . 
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With  these  values,  the  power  supplied  to  the  antenna  per  unit  current  is 

ZT  /|l(0)]2  =  R  (k)  =  1.71  (1  +  140  a jf> ) watt s / amp 2 ;  (116a) 

the  power  transferred  beyond  the  radian  sphere  per  unit  current  is 

2Tftr/jl(0)|^  =  1.71(1  f  2a/p)  watts/amp2  ;  (116b) 

and  the  power  dissipated  as  heat  in  the  medium  within  the  radian  sphere  is 

=  236  a/p  watts/amp*"  .  (116c) 

When  a  =  P  ,  the  corresponding  values  are 


o  Vr*  / 1  •>  i  a  i  1  ^  *>  m  %  i  a  r  d  t  /  2 

£le/|i\u;j  -  n.  (kj  =  y  wdub/diuu 

/  1  1  V 

Vina; 

_  2  2 

?T  /Irm  1  —  0  8  1  WO  f  f  o  /awn 

~-er'  I-'"'*  - , 

(117b) 

2 

P  =  1989-7  watts/amp 
r 

(117c) 

From  these  numerical  results  it  is  clear  Shat  when  a/p  is  as  small  as 
10  or  smaller,  90  percent  or  more  of  the  power  supplied  to  the  antenna  is 
transferred  beyond  the  radian  sphere.  When  ci/p  is  no  greater  than  0.  1,  only 
about  o  percent  is  dissipated  outside;  92  percent  is  used  to  heat  the  medium 
inside  the  radian  spnere.  When  a  =  P  virtually  all  of  the  power  is  dissipated 
as  heat  within  the  radian  sphere.  The  fraction  transferred  beyond  it  is  only 
about  0.  04  percent.  Note  that  when  a  =  P  the  input  susceptance  of  the  antenna 
is  small,  so  that  virtually  the  entire  impedance  is  resistive.  The  antenna 
effectively  does  not  radiate,  but  acts  like  a  pair  of  electrodes  with  very  small 
surface  area- -hence  the  rather  high  resistance. 

Conclusion 

The  general  problem  of  a  cylindrical  antenna  immersed  in  a  dissipative 
medium  has  been  formulated  in  a  manner  that  permits  the  determination  not  only 
of  the  distribution  of  current  and  the  admittance,  but  also  of  the  electromagnetic 
field.  The  analytical  procedure  is  approximate  but  quantitatively  sufficiently 
accurate  to  be  of  practical  value.  It  may  be  extended  to  treat  coupled  antennas 
in  a  dissipative  medium. 
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